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We report on simulations of mildly turbulent convection in spherical wedge geometry with varying density stratification. 
We vary the density contrast within the convection zone by a factor of 20 and study the influence of rotation on the 
solutions. We demonstrate that the size of convective cells decreases and the anisotropy of turbulence increases as the 
stratification is increased. Differential rotation is found to change from anti-solar (slow equator) to solar-like (fast equator) 
at roughly the same Coriolis number for all stratifications. The largest stratification runs, however, are sensitive to changes 
of the Reynolds number. Evidence for a near-surface shear layer is found in runs with strong stratification and large 
Reynolds numbers. 
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1 Introduction 

Numerical simulations of turbulent convection in spherical 
geometry have become a standard tool in the study of differ- 
ential rotation and magnetism in the solar and stellar context 
(Miesch & Toomre 2009). The current state of the art mod- 
els can reproduce many aspects of the solar internal rota- 
tion (e.g. Miesch et al. 2006), and large-scale oscillatory dy- 
namo action occurs when rotation is rapid enough (Brown 
et al. 2010, 201 1). However, reproducing the solar cycle has 
turned out to be elusive, even though large-scale oscillatory 
fields are now seen in some simulations with the solar rota- 
tion rate, too (e.g. Ghizaru et al. 2010; Racine et al. 2011). 
The reason for the remaining discrepancies may lie in the 
fact that much of the physics have either to be simplified 
or neglected altogether due to severe numerical constraints 
(e.g. Kapyla 2011). Furthermore, the simulations that are 
being carried out are so demanding that often only a single 
or a few representative cases can be done. Although many 
results, such as the change from anti-solar (slow equator) to 
solar-like (fast equator) differential rotation as the rotation 
rate increases (e.g. Chan 2010; Kapyla et al. 2011a), and 
the appearance of mostly axisymmetric large-scale mag- 
netic fields (e.g. Oilman 1983; Glatzmaier 1985, Browning 
et al. 2006; Brown et al. 2010, 201 1; Kapyla et al. 2010) ap- 
pear robust, their exact dependence on different simulation 
parameters has not been explored in detail. 

Here we study the effect of density stratification on ro- 
tating spherical shell convection. Our main goal is to study 
how the transition from anti-solar to solar-like rotation is af- 
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fected. This is relevant because most convection models in 
spherical shells still have rather modest density stratification 
in comparison to the Sun. Although most of the mass within 
the convection zone is located near the base, fast downflows 
at the vertices of convection cells originate near the surface. 
The effect of these downflows on angular momentum trans- 
port is yet unclear. We are also interested in the statistical 
properties, such as anisotropy, of turbulence as the stratifi- 
cation is increased. 

2 Model 

Our model is based on that used by Kapyla et al. (2010, 
2011a). We model a segment of a star, i.e. a "wedge", in 
spherical polar coordinates, where (r, 9, (p) denote the ra- 
dius, colatitude, and longitude. The radial, latitudinal, and 
longitudinal extents of the computational domain are given 
by 0.7R < r < R, 9o < 9 < TT - do, and < 4) < (f>Q, 
respectively, where R is the radius of the star. In all of our 
runs we take 9o — tt/S and (j>o = 7t/2. 

We solve the following equations of compressible hy- 
drodynamics in a frame of reference rotating with angular 
velocity 

Dlnp 

Dt 
Du 

lyt 

Ds 



= V • u, 



1 



= g 2r2o X M + - ( V • 2iypS - Vp) , 



(1) 



(2) 



(3) 



where D / Dt = d / dt + u ■ 'V is the advective time deriva- 
tive, p is the density, u is the velocity, s is the specific 
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entropy, T is the temperature, and p is the pressure. The 
fluid obeys the ideal gas law with p — (7 — l)pe, where 
7 — cp/cv — 5/3 is the ratio of specific heats at constant 
pressure and volume, respectively, and e = cyT is the in- 
ternal energy. 

Furthermore, ly is the kinematic viscosity, K is the ra- 
diative heat conductivity, xt is the unresolved turbulent heat 
conductivity, and g is the gravitational acceleration given by 

9 = — ^I^-*^' ("^^ 

where G is the gravitational constant, M is the mass of the 
star, and f is the unit vector in the radial direction. We omit 
the centrifugal force in our models. The rate of strain tensor 
S is given by 

S^j = ^{Ui-j + Uj-i) - V ■ M, (5) 

where the semicolons denote covariant differentiation; see 
Mitra et al. (2009) for details. Unlike in our previous studies 
(Kapyla et al. 2010, 20 11 a), we omit stably stratified layers 
below and above the convectively unstable layer. 

2.1 Initial and boundary conditions 

In the initial state the atmosphere is adiabatic and the hy- 
drostatic temperature gradient is given by 

— = -g (6) 

dr cv(7 — l)(m + 1) ' 

where m = 1.5 is the poly tropic index. We use Eq. ^ as the 
lower boundary condition for the temperature. This gives 
the logarithmic temperature gradient V (not to be confused 
with the operator V) as 

V = 51nr/91nj3= (m + l)~^ (7) 

Density stratification is obtained by requiring hydrostatic 
equilibrium. The heat conduction profile is chosen so that 
radiative diffusion is responsible for supplying the energy 
flux in the system, and K decreases rapidly within the con- 
vection zone (see, Fig.[Tll. 

The radial and latitudinal boundaries are taken to be im- 
penetrable and stress free, according to 
du0 _ ue_ du^ _ U0 
dr r ' dr r 



Ur = 0, 



dUr 

'do 



= = 0, 



= ua, cot ( 



{r = 0.7R,R), (8) 
9 = 9o,7r-eo). (9) 



On the latitudinal boundaries we assume that the thermody- 
namic quantities have zero first derivative, thus suppressing 
heat fluxes through the boundaries. 

On the upper boundary we apply a black body condition 



or 



(10) 



where a is the Stefan-Boltzmann constant. In our runs we 
use a modified value for cr that takes into account that our 
Reynolds and Rayleigh numbers are much smaller than in 
reality, so K is much larger and therefore the flux too high. 
The black body boundary for the temperature has previ- 
ously been used in mean-field models of Brandenburg et al. 
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Fig. 1 Profiles of density (top panel), number of pressure 
scale heights, iVp — \n ^, where po is the pressure at r = 
0.7R (middle), and xt and K (bottom). In the bottom panel 
1/ and p from Set A are used as normalization factors. 



(1992). In our runs K is negligibly small near the surface so 
that the unresolved convective energy flux transports prac- 
tically all of the energy through the upper boundary. This is 
similar to what is commonly used in the ASH simulations 
(e.g. Brun et al. 2004). 

2.2 Dimensionless parameters 

We obtain non-dimensional quantities by choosing 

R = GM = po = cp = 1 , (11) 

where po is the density at 0.7 R. The units of length, veloc- 
ity, density, and entropy are then given by 



[x]=R, [u] = ^GM/R , [p]=po, [s] = cp. (12) 

The simulations are governed by the Prandtl, Reynolds, 
Coriolis, and Rayleigh numbers, defined by 

2O0 



Pr = 



Ra 



Re = 



Xt vkf 
GM{Ar)^ / 1 ds 
i^XtR'^ \ cp dr 



Co = 



^rms^'f 



(13) 



(14) 



where xt is the turbulent thermal conductivity in the middle 
of the convection zone (i.e. at r^ — 0.85i?), fcf = 27r/Ar 
is an estimate of the wavenumber of the energy-carrying 
eddies, Ar = 0.3R is the thickness of the layer, and 

liims — ^ + ^e) is the rms velocity, where the an- 
gle brackets denote volume averaging. In our definition of 
Urms we omit the contribution from the 0-velocity, because 
its value is dominated by effects from the differential ro- 
tation (Kapyla et al. 2011a). Sometimes we show Wrms('") 
which is the fluctuating rms velocity as a function of radius 
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Fig. 2 Flux balance from Run B3. The different contribu- 
tions are due to radiative diffusion (solid line), resolved con- 
vection (dashed), unresolved turbulence (dotted), flux of ki- 
netic energy (dot-dashed), and viscosity (triple-dot-dashed). 
The red dotted lines denote the zero level and the total lu- 
minosity through the lower boundary. 



and from which we have subtracted the azimuthally aver- 
aged velocities. The entropy gradient, measured at r„^, is 
given by 



_ J_ds\ 
cp dry, 



V„, - V 



ad 



Hp 



(15) 



where Vm = {dlnT/d\np)r^, and Hp is the pressure 
scale height at r,n. Due to the fact that the initial stratifica- 
tion is isentropic, we quote values of Ra from the thermally 
saturated state of the runs. 

The energy that is deposited into the domain at the base 
is controlled by the luminosity parameter 

Lo 



P0(GM)3/2i?l/2 



where Lq = AnrfFi, is the constant luminosity, and F^, — 
— {KdT/dr)\r=o.jR is the energy flux imposed at the bot- 
tom boundary. Furthermore, the stratification is determined 
by the normalized pressure scale height at the surface 



(7 - ^)cvTi 



(17) 



GM/R 

where Ti — T{r — R). Similar parameter definitions were 
used by Dobler et al. (2006). We use three different values 
(0.09, 0.02, 8 • 10^3) of ^ which result in density conti-asts 
of 5, 30, and 10^, respectively (see Fig.[T]). Now the convec- 
tion zones span between roughly 2.5 and 7.5 pressure scale 
heights. 

The simulations were performed using the PENCIL 
CODllil which uses sixth-order explicit finite differences 
in space and a third-order accurate time stepping method; 
see Mitra et al. (2009) for further information regarding the 
adaptation of the PENCIL Code to spherical coordinates. 



3 Results 

We have performed four sets of simulations differing by 
their density stratification (see Table [U and Reynolds num- 
ber. We vary the rotation rate within each set so that the 
Coriolis number changes by roughly an order of magnitude. 
We increase the gravity by a factor of in Sets C and D in 
order to limit the Mach number to roughly 0. 1 near the sur- 
face. The grid resolution in Sets A-C is 128 x 256 x 128. In 
Set C this means that the ratio of the pressure scale height to 
the radial grid spacing at the surface is Hp /5r m 3.4, which 
is on the Umit of resolving the structure. We have remeshed 
snapshots from the saturated states of the runs in Set C to 
double resolution (Set D) where we are also able to increase 
the Reynolds number. 

3.1 Flux balance 

In contrast to our earlier studies using a polytropic setup 
with m — \ (Kapyla et al. 2010, 2011a), we now use a 
setup in which convection transports the majority of the 
flux. This is achieved by decreasing the heat conductivity K 
within the convection zone and introducing a turbulent heat 
conductivity xt which is responsible for unresolved convec- 
tive transport of heat (e.g. Chan & Sofia 1996; Brun et al. 
2004). We apply a constant value of Xt, of the order of the 
kinematic viscosity v, in the bulk of the convection zone 
(G.75i? < r < 0.98i?) and an order of magnitude larger 
value above r > 0.98R in order to transport the flux through 
the upper boundary. Below r — 0.75R, Xt goes smoothly to 
zero, see Fig.[T] 

To verify that the system is in thermal equilibrium, we 
consider the radiative, convective, kinetic, viscous, and tur- 
bulent energy fluxes, defined as 



(16) Fr 



ad 



dr ' 



-cppu'^T', 



Fx 



kin 



F^ 



F 



turb 



(18) 

(19) 
(20) 

(21) 
(22) 



http : / /pencil -code . googlecode . com/ 



_- ds 

--PTX.Q-^. 

where the averages are taken over and 0. Representative 
results from Run B3 are shown in Fig.|2l Radiative diffusion 
transports the total flux through the lower boundary and de- 
creases rapidly as a function of r. The radiative flux is less 
than 10 per cent above r = 0.85i?. The flux due to resolved 
convection is responsible for transporting the majority of 
the luminosity within the convection zone. The flux of ki- 
netic energy is directed downwards and is responsible for 
roughly 10 per cent of the flux near the surface. Note that 
the maxima of Fconv and Fy\n are significantly larger in the 
non-rotating cases. The unresolved turbulent flux is small 
in the bulk of the convection zone and carries the flux out 
through the outer boundary. The viscous flux is small in all 
of our runs. The flux balance is similar to the ASH sim- 
ulations (e.g. Brun et al. 2004; Miesch et al. 2008) which 
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Table 1 Summary of the runs. Here, Ma = Urms/ \J GM/R, fcjj — Ail/f^cq, A51 — 51cq — f^poio, where ficq = 
^{R^ = 7r/2) and ftpo\c = ^{R, & — Sq). Furthermore, Ey^ — {^pu^) is the volume averaged total kinetic energy, and 
Em — ^(/o(?I^ + Ug)) and = ■^{pu'^) are the kinetic energies of the meridional circulation and differential rotation, 
respectively. The coefficients lui, ws, and W5 represent the expansion coefficients in Eq. ( l26b and will be discussed in 
Sect.[33]below. 
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Fig. 3 Radial velocity Ur, normalized by the local sound speed from r = 0.98i? for runs in Sets A (top row), B, C, 
and D (bottom row). The rotation rate increases from left to right. The longitude extent has been duplicated fourfold for 
visualization purposes. 
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Fig. 4 Radial velocity Ur in the meridional plane (p ~ (fiQ 
from Runs AO, BO, and DO in units of ^JGM/R. 

employ stratification from a ID solar model and somewhat 
different profiles of the diffusion coefficients. 

3.2 Properties of convection 

VisuaUzations of the radial velocity near the surface of the 
star for all of our models are shown in Fig. [3] We find that 
as the stratification increases the size of convection cells de- 
creases. This is a consequence of the decreasing pressure 
scale height near the surface. In the case of the highest strat- 
ification. Sets C and D, the granulation pattern is similar 
to the high-resolution run reported by Miesch et al. (2008) 
with a comparable stratification. As noted above, the resolu- 
tion in Set C is close to critical in resolving the stratification 
properly which is also manifested by numerical artefacts in 
Fig. [3] The higher resolution runs in Set D, however, are 
well behaved and show similar convection patterns in the 
non-rotating and slowly rotating cases. 

Furthermore, as the rotation rate is increased, one sees 
the formation of a cartridge belt-like pattern that is also 
known as banana cells. However, these structures become 
less pronounced at larger stratification in Set C. In Set D, on 
the other hand, strong banana cells are again observed. The 
main difference between Sets C and D is that in the latter the 
Reynolds number is modestly increased. The smaller size 
of convection cells at high stratification leads to a smaller 
"effective" Reynolds number based on the horizontal eddy 
size, which is not well reflected by our definition of Rc. It 
is possible that in Run C4 this effective Rc is below critical 
to excite the formation of banana cells and strong prograde 
differential rotation. 

Another way to see the difference between weak and 
strong stratification is to visualize the flows in the merid- 
ional plane. In Fig.|4]we show a cut of the radial velocity at 
= 00 from Runs AO, BO, and DO in the saturated state. 
Whereas the downflows in Run AO go easily through the 
whole convection zone, smaller-scale structures originating 
from the surface appear already in Run BO. For the strongest 
stratification the anisotropy of the flow is clear also to the 
naked eye with smooth large-scale flows in the deep lay- 
ers and small-scale irregular flows near the surface. How- 
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Fig. 5 From top to bottom: the three uppermost panels 
show the radial dependence of the fluctuating velocity com- 
ponents averaged over Q and for nonrotating Runs AO, 
BO, and CO, respectively. The fourth and fifth panel show 
the turnover time Tto = -ffp/urms('') divided by tq = 
(uimsofcf )"^. and the vertical anisotropy parameter Ay for 
the same runs as indicated by the legends. Here Urmso = 

\^ \{u1 + Ug) is used as a normalization factor. 

ever, the strongest downflows are still able to go all the way 
through the convection zone. 

The squares of the fluctuating velocity components are 
shown in Fig. |5] For the weakest stratification the profiles 
are almost symmetrical with respect to r,„. The apparent 
asymmetry of the Q and velocities is likely due to the 
large size of the convective cells in comparison to the do- 
main size. For larger stratification the velocities increase 
near the surface and the anisotropy of the horizontal veloc- 
ities decreases. An important effect that follows is that the 
convective turnover time, defined as 

Tto = i?p/Urms(f), (23) 

where i/p is the local pressure scale height, changes sub- 
stantially between the bottom and the surface. In the Run AO 
Tto is almost constant in the whole layer whereas in Runs BO 
and CO it varies by factors of 9.1 and 31, respectively. Pro- 
vided that mixing length arguments hold, the rotational in- 
fluence on the flow, measured by the Coriolis number, has 
the same variation as a function of radius. 
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A. 



We define the vertical anisotropy parameter as 



(24) 



where the averages are taken over 9 and (j), and the primes 
denote that ^-averaged mean velocities are subtracted. For 
weak stratification, Ay is at most 0.4 in the middle of the 
convection zone. Note that Ay = —1 at the boundaries 
due to the impenetrable boundary conditions. For Runs BO 
and CO, the turbulence is more anisotropic and the point at 
which Ay changes from negative to positive moves closer 
to the boundaries. The anisotropy measured by Ay also in- 
creases and in Runs BO and CO it is more than double the 
value of Run AO. In comparison to forced turbulence sim- 
ulations of Brandenburg et al. (2011a), our results corre- 
spond best to cases with small scale separation, i.e. large- 
scale forcing. This is consistent with large convective cells 
that span the whole depth of the convection zone. The im- 
portance of Ay is that in rotating turbulence, it acts as a 
source for the vertical A-effect (e.g. Riidiger 1989; Kapyla 
& Brandenburg 2008) which drives radial differential rota- 
tion. 

3.3 Differential rotation 

The differential rotation profiles from all runs with flo ^ 
are shown in Fig. |6] We find that an anti-solar differential 
rotation pattern with strong meridional circulation forms at 
the lowest rotation rates. As increases, equatorial accel- 
eration gradually develops. However, in many cases (e.g. 
Runs B3, B4, C4, D4) there is a minimum of the local an- 
gular velocity, fl — u^/r sin 9, at mid-latitudes and a polar 
vortex at high latitudes. Similar profiles have been reported 
by Miesch et al. (2000) and ElHott et al. (2000). Large-scale 
vortices arise in Cartesian convection simulations at suffi- 
ciently high rotation rates (Chan 2007; Kapyla et al. 2011b; 
Mantere et al. 2011). It is unclear whether the polar vortices 
in spherical geometry are related to the vortex instability but 
it is an intriguing possibility. 

We quantify the horizontal differential rotation by the 
parameter 



(25) 



where flp = ^[n{R,9o) + n{R,TT ~ 9q% and f^oq = 
ri(i?, 7r/2). The results for Sets A to D, along with cor- 
responding data from Kapyla et al. (201 la) are shown in 
Fig. [8] We find that for Co « 1, fcn is the largest for the 
smallest stratification. The results for Sets B and D seem 
to converge for more rapid rotation and produce solar-like 
rotation (fc^ > 0) for Co > 5. 

In Set C, however, the transition to a solar-like profile 
does not yet occur in the parameter range studied here, al- 
though the largest Coriolis number is of the order of 7. An- 
other factor that comes into play is the fact that the effective 
Reynolds number based on the typical scale of convection 
cells is reduced in the runs with the largest stratification. 
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Fig. 6 Rotation profiles = u^j (r Wi9) + 57o, nor- 
malised by f^o, from all runs with ilo 7^ 0. 



This seems to be confirmed by the simulations in Set D, 
which are the higher Reynolds number counterparts of the 
runs in Set C, although for Run D4, fco is still negative. This 
is surprising given the profile seen in Fig. |6] which clearly 
shows a rapidly rotating equator. The discrepancy is due to 
a sharp negative radial gradient of near the surface at the 
equatorial regions in Run D4 (see Fig. |7]i. This is similar 
to the near-surface shear layer observed in the Sun (e.g., 
Benevolenskaya et al. 1999). 

The formation of a negative near-surface shear layer 
only occurs for strong stratification. Furthermore, Runs C4 
and D4 (where Re is twice as large), suggests that negative 
near-surface shear also requires large Reynolds numbers. It 
is still unclear whether the current simulations can really 
capture the physics of the solar near-surface shear layer (e.g. 
Miesch & Hindman 201 1), but the present results might in- 
dicate a path that is worth following. Measuring fcji from a 
little deeper down at r = 0.95i? gives 0.06 which is similar 
to the results from Runs B3 and B4 (see TablelTJ. 

We also note that the results from Kapyla et al. (201 la) 
roughly fall in line with Sets B and D for Co < 3. However, 
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Fig. 7 Rotation profiles at the equator from Runs A4 
(solid line), B4 (dashed), C4 (dot-dashed), and D4 (triple- 
dot-dashed). 
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Fig. 8 Differential rotation parameter for Sets 
as indicated by the legend. The black triangles refer 
taken from Kapyla et al. (20 11 a). 
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in the rapid rotation regime, /co from Kapyla et al. (20 11 a) 
is consistently larger than in the current results with the ex- 
ception of Set A. This is probably due to the difference in 
the setups, i.e. we omit here the stably stratified overshoot 
layer below the convection zone and the isothermal cooUng 
layer near the surface. 

It is customary to represent the latitudinal profile of the 
angular velocity in terms of Gegenbauer polynomials, i.e., 

ijJiPI (cos &) I sin 9. 



E 

e=i,3,5,... 



(26) 



(We use here a definition where all associated Legen- 
dre polynomials with odd m values are positive, i.e.. 
Pi (cos 6) = sin 0, for example.) We have determined the 
coefficients wi, ^3, and 015 via a fitting procedure using 
sin^6' as weighting factor to put more emphasis on the equa- 
torial regions. The resulting coefficients are given in the last 
3 columns of Table[T] Note that ^3 changes sign for Co > 3. 
This is also consistent with Fig. |9] where we plot the lati- 
tudinal profiles of 2il/urnisk[ together with their fits from 
Eq. ( l26b . The change of sign of ut^ gives a more robust in- 
dicator of the transition from antisolar to solar-like rotation 
than just the fco parameter. 
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Fig. 9 Fits to the latitudinal differential rotation at r = 
0.98i? for Sets A (soHd), B (dotted), C (dashed), and D 
(dash-dotted). The actual data are shown as thin (red) lines. 

4 Conclusions 

We study turbulent convection in spherical shells with vary- 
ing density stratification. We find that the typical size of 
convection cells decreases as the stratification increases, 
which is in accordance with mixing length arguments. At 
the same time the anisotropy of turbulence increases and 
the turnover time varies by more than an order of magni- 
tude from the base to the top of the convection zone. 

Although convection seemingly changes greatly as the 
stratification increases, the rotation profiles and their qual- 
itative trend as a function of the Coriolis number change 
surprisingly little. However, we find that the results for the 
largest density stratification are sensitive to changes of the 
Reynolds number and that apparently smaller latitudinal 
differential rotation for large Coriolis numbers is obtained 
for large stratification. Measuring the differential rotation 
simply as a difference between the surface values of il at 
the equator and at high latitudes turns out to give mislead- 
ing results for our high-resolution runs with the largest strat- 
ification. This is due to the self-consistent generation of a 
sharp radial gradient of fl near the surface, reminiscent of 
the near-surface shear layer in the Sun. A similar feature is 
discernible in the highly stratified simulations of Bessolaz 
& Brun (20 11). However, it remains to be seen whether this 
is a robust feature. 

Another interesting aspect of increasing stratification is 
related to the generation of large-scale magnetic fields, be- 
cause some of the contributions to the a-effect of mean-field 
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dynamo theory are proportional to density stratification (e.g. 
Krause & Radler 1980). Furthermore, the negative magnetic 
pressure instability (e.g. Brandenburg et al. 201 lb; Kapyla 
et al. 201 Ic), that can lead to magnetic field concentrations 
of the form of active regions, becomes stronger when strat- 
ification increases. We plan to revisit these issues in forth- 
coming papers. 

Acknowledgements. We thank the referee for the suggestion to 
add the expansion in terms of Gegenbauer polynomials. Compu- 
tational resources granted by CSC - IT Center for Science, who 
are financed by the Ministry of Education, and financial support 
from the Academy of Finland grants No. 136189, 140970 (PJK), 
218159 and 141017 (MJM), and the European Research Council 
under the AstroDyn Research Project 227952 (AB) are acknowl- 
edged. 



Mitra, D., Tavakol, R., Brandenburg, A., Moss, D.: 2009, ApJ 697, 
923 

Racine, E., Charbonneau, P., Ghizaru, M., Bouchat, A., Smo- 

larkiewicz, P K.: 201 1, ApJ 735, 46 
Rudiger, G.: 1989, Dijferential Rotation and Stellar Convection: 

Sun and Solar-type Stars (Akademie Verlag, Berlin) 



References 

Benevolenskaya, E.E., Hoeksema, J.T., Kosovichev, A.G., Scher- 
rer, PH.: 1999, ApJ 517, L163 

Bessolaz, N., Brun, A.S.: 2011, ApJ 728, 115 

Brandenburg, A., Moss, D., Tuominen, I.: 1992, A&A 265, 328 

Brandenburg, A., Radler, K.-H., Kemel, K.: 2011a, A&A (submit- 
ted), arXiv: 1108.2264 

Brandenburg, A., Kemel, K., Kleeorin, N., Mitra, D., Ro- 
gachevskii, 1.: 2011b, ApJL 740, L40 

Brown, B.P, Browning, M.K., Miesch, M.S., Brun, A.S., Toomre, 
J.: 2010, ApJ 711, 424 

Brown, B.P, Miesch, M.S., Browning, M.K., Brun, A.S., Toomre, 
J.: 2011, ApJ 731, 69 

Browning, M.K., Miesch, M.S., Brun, A.S., Toomre, J.: 2006, 
ApJ 648, L157 

Brun, A.S., Miesch, M.S., Toomre, J.: 2004, ApJ 614, 1073 

Chan, K.-L., Sofia, S.: 1996, ApJ 466, 372 

Chan, K.-L.: 2007, AN 328, 1059 

Chan, K.-L.: 2010, lAUS 264, 219 

Dobler, W., Stix, M., Brandenburg, A.: 2006, ApJ 638, 336 
Elliott, J.R., Miesch, M.S., Toomre, J.: 2000, ApJ 533, 546 
Ghizaru, M., Charbonneau, P., Smolarkiewicz, P.K.: 2010, ApJ 

715,L133 
Gilman, PA.: 1983, ApJS 53, 243 
Glatzmaier, G.A.: 1985, ApJ 291, 300 
Kapyla, PJ., Brandenburg, A.: 2008, A&A 488, 9 
Kapyla, P.J., Korpi, M.J., Brandenburg, A., Mitra, D., Tavakol, R.: 

2010, AN 331, 73 
Kapyla, P.J., Mantere, M.J., Guerrero, G., Brandenburg, A., Chat- 

terjee, R: 2011a, A&A, 531, A162 
Kapyla, PJ., Mantere, M.J., Hackman, T: 201 lb, ApJ , 742, 34 
Kapyla, P.J., Brandenburg, A., Kleeorin, N., Mantere, M. J., Ro- 

gachevskii, 1.: 2011c, MNRAS (submitted), arXiv: 1104.4541 
Kapyla, PJ.: 201 1, AN, 332, 43 

Krause, R, Radler, K.-H.: 1980, Mean-field Magnetohydrodynam- 

ics and Dynamo Theory, Pergamon Press, Oxford 
Mantere, M.J., Kapyla, P.J., Hackman, T: 2011, AN (submitted), 

arXiv: 1109.43 17 
Miesch, M.S., Elliott, J.R., Toomre, J., et al.: 2000, ApJ 532, 593 
Miesch, M.S., Brun, A.S., Toomre, J.: 2006, ApJ 641, 618 
Miesch, M.S., Brun, A.S., DeRosa, M.L., Toomre, J.: 2008, 

ApJ 673, 557 
Miesch, M.S., Toomre, J.: 2009, AnRFM, 41, 317 
Miesch, M.S., Hindman, B.W.: 2011, ApJ (submitted), 

[arXiv: 1106:41071 



© 2011 WILEY- VCH Verlag GmbH &Co.KGaA,Weinheim 



www.aii-joumal.org 



